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1 Abstract 

> ; 

We report on the occurrence of negative-phase velocity (NPV) planewave propagation 
(y-j \ in the ergosphere of a rotating black hole. By implementing the Kerr metric, it is 

demonstrated that regions of NPV propagation are concentrated at the equator of 
the ergosphere, while NPV propagation is less common towards the polar regions. 
', Increasing the angular velocity of the black hole exaggerates the NPV concentration 

"^h ' at the equator. NPV propagation is not observed outside the stationary limit surface. 

a.' 
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1 Introduction 

The term negative-phase velocity (NPV) refers to a mode of electromagnetic planewave 
propagation in which the phase velocity vector is directed opposite to the direction of time- 
averaged energy flux [1, 2]. Of the many interesting consequences of NPV propagation, it is 
the phenomenon of negative refraction which has recently attracted much attention in the 
electromagnetics/optics and materials research communities [3]. This interest may be traced 
back to the first experimental reports of negative refraction in 2001, involving the microwave 
illumination of an artificial metamaterial [4] . 

We have recently begun to examine the possibility of gravitationally-assisted NPV prop- 
agation in vacuum [5, 6]. In the absence of a 'material', the source of the NPV may be 
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provided by the curvature of spacetime. The electrodynamic formulation adopted is founded 
on the formal analogy between (a) electromagnetic propagation in curved spacetime in free 
space; and (b) electromagnetic propagation in flat spacetime in a fictitious, instantaneously 
responding bianisotropic medium. This approach — which was originally put forward by 
Tamm and subsequently described by numerous authors [7]— [13] - - enables wave propaga- 
tion in curved spacetime to be examined using standard electromagnetic techiques. Thereby, 
the feasibility of NPV propagation for a general gravitomagnetic metric has been established 
[14]- 

We further develop the concept of gravitationally-assisted NPV propagation in this paper 
with a detailed numerico-theoretical study based on a particular gravitomagnetic metric, 
namely the Kerr metric [15]. Specifically, we explore the spacetime surrounding a rotating 
black hole, in search of regions where NPV propagation occurs. Our search is restricted to the 
physically probeable spacetime outside the outer event horizon. The theoretical framework 
is described in Section 2, numerical results are presented and analyzed in Section 3, and 
conclusions are itemized and discussed in Section 4. 



2 Theoretical framework 

The electrodynamics of the classical vacuum in a curved spacetime provides the setting 
for our numerico-theoretical study. Thus, we consider the covariant and the contravariant 
electromagnetic field tensors f a p and h a ^, respectively, which appear in the source-free 
covariant Maxwell equations^ 



fa/3;u + f/3i>;a + fva-fi ~ 



h 



a/3 



;/3 







where the subscript - u indicates the covariant derivative with respect to the uth space- 
time coordinate x v . The spacetime curvature is described by the metric g a p with signature 

(+,-,-,-)■ 



2.1 Electrodynamic formulation 

We follow the standard approach [8]— [11] - originally proposed by Tamm 
covariant equations (1) are replaced by their noncovariant equivalents 



fa/3,u + f(3v,a + fva,P 



o 



o 



wherein the 



(2) 



§ Greek indexes take the values 0, 1, 2 and 3; Roman indexes take the values 1, 2 and 3; x° — ct where c 
is the speed of light in vacuum in the absence of all gravitational fields; whereas x 1 ' 2,3 are the three spatial 
coordinates. 



2 



with g = det [g a p\ and the subscript >v denoting ordinary differentiation with respect to the 
z/th spacetime coordinate. More conveniently, we express (2) in the familiar form 



m,n 

De,e = , — Dtp + Et mn H m ^ 



(3) 



(4) 



where e^ mn is the three-dimensional Levi-Civita symbol. The conventional electromagnetic 
field vectors E_ , B_ , D and H_ are the 3 vector-equivalents associated with the respective 
components 

Ee = feo , 

Be = (l/2)e£ mn f rnn 
D e =(-gf 2 h*\ 
H e =(l/2)e emn (-g) 1/2 h mn 

introduced in (3). Thus, in (3), B^g and Dgg represent the divergences of B_ and D_, whereas 
the £th vector component of the curls of E_ and H_ are represented as Ee mn E m)n and Eg mn H m ^ ni 
with the derivatives with respect to time of the £th vector component of B_ and D_ being repre- 
sented as B/ifl and D^ . The electromagnetic field vector components satisfy the constitutive 
relations 

De = ~it m E m + Ei mn T m H n 

Bg / j£ m H rn £lmn Tm E n 

where the tensor 7£ m components and vector T m components are defined as 

li m = ~ {-9) , 

9oo } . (6) 

r jy). 



(5) 



flOm 



#00 



Thus, in terms of the noncovariant formulation (2), vacuum is represented by the fictitious, 
instantaneously-responding, bianisotropic medium characterized by (5). 
For later use, we recast (3) and (5) in 3 x 3 dyadic/ 3 vector form as 



and 



d 

V x E(ct,r) + -B(d,r) = 

ot 

d 

V x H(ct,r) - -^D(ct,r) = 



D(ct,r) = e ( ,7(ct,r) • E(ct,r) - -£(ct,r) x H_(ct,r) 
= c 

B(ct,r) = Hol(ct,r) • Hict,r) + -£(ct,r) x E(ct,r) 
= c 



(7) 



(8) 



respectively, wherein 7(ct,r) is the dyadic-equivalent of 7^ m , T(ct,r) is the vector-equivalent 
of r m ; the scalar constants e and /i denote the permittivity and permeability of vacuum in 
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the absence of a gravitational field; and SI units are adopted. In moving from (1) to (7) and 
(8), the spacetime coordinate x a has been separated into space r and time t coordinates. 
The vector cross products and dyadic-vector dot products in (7) and (8) apply to vectors 
and dyadics associated with the three spatial dimensions of flat spacetime [16]. However, it 
is emphasized that the spacetime underlying (7) and (8) remains curved. 



2.2 Kerr spacetime 

We consider a specific gravitomagnetic spacetime, namely the Kerr spacetime [15], which 
describes a rotating black hole. In Cartesian coordinates — r = (x, y, z) — the components 
of the Kerr metric are given as [17] 



#oo 

#01 

#02 
#03 

#11 
#12 
#13 
#22 
#23 



1-A, 
. Rx 



= -A 



R 2 + «L 



R 2 + a 2 



bh 



= -A 



1-A 



Rx + a lbh y 



-A 



-A 



-1-A 



R 2 + «L 

(Rx + a Ihh y) (Ry - a rbb x) 

(R 2 + <J 2 

(Rx + a Ihh y) z 

(R 2 + a?J R ' 
a Thh x - Ry 



(9) 
(10) 

(11) 
(12) 

(13) 
(14) 
(15) 
(16) 
(17) 
(18) 

(19) 

" ~ i? 4 + (a lbh zf ' (20) 

The components (9)— (18) of the metric describe a black hole of geometric mass m rbh , rotating 
about the Cartesian z axis, in conventional units. The term a Ihh is a measure of the black 
hole's angular velocity (sometimes referred to as the specific angular momentum). 



R 2 



= -A 



(Ry - a rbh x) z 

(R 2 + «L) R ' 
#33 = -i-A (|) 2 , 

with g a p = gp a . The quantity R is given implicitly via 



R 2 



also, 



2,2,2 2 

x +y +z -a rbh 



2m Tbh R 3 
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We note that the Kerr metric is often represented in terms of Boyer-Lindquist coordinates 
[17, Sec. 19.3]. However, for the planewave analysis undertaken in the following sections, 
the Cartesian representation (9)-(18) is more convenient. 

Using the definitions (6) together with (9)— (18), we obtain 



7n 

712 
713 
722 

723 
733 



6 [R 4 + (a lhh z) 2 ] 



1 - A 



Rx + a Ibh y 
R 2 + a 2 



bh 



3 (Rx + a lhh y) (Ry - a lbh x) 
— Zo m rhh ri — 



(R 2 + <J 



-26 m Ihh R z—— : — - — 



R 2 



6 [R 4 + (a Ihh z) 2 } 



1 - A 



a ibh x Ry 

R 2 + a 2 



or D 2 Ril a rbh x 

-25m^R z R2 + 



bh 



6 [R 4 + (a Ihh z) 2 ] 



1 - A 



R 



with 7^ m = 7 m£ , and 



-26m Ibh R 
—25m rhh R 



3 i Rx + a rbh y 

n 



R 2 + < h 



3 i Ry ci r bh x 



wherein 



T 3 = -26 m Ibh R z , 
6 = 



(21) 
(22) 
(23) 

(24) 

(25) 
(26) 



(27) 

(28) 
(29) 

(30) 

R 4 + (a lhh z) - 2m lhh R 3 

We concentrate on the regime m 2 hh > a 2 hh . Furthermore, as our interest lies in physically 
observable phenomenons, we restrict attention to the region R > R + where 



R+ = m lbh + a/tti^ 



(31) 



represents the outer event horizon of the black hole. Finally, we note that R + < Rs + where 



Rs+ = m lhh + \ m 2 bh - 



(l rh uZ 



R 



s+ 



(32) 



delineates the stationary limit surface of the black hole. 
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2.3 Planewave propagation 

Let us consider the propagation of plane waves in the medium characterized by the con- 
stitutive relations (8). Our strategy is to focus upon an arbitrary neighbourhood 1Z whose 
location is specified by the Cartesian coordinates (x,y,z). The neighbourhood is presumed 
to be sufficiently small that the nonuniform metric g a p may be approximated by the uniform 

metric g a p throughout 1Z. Thus, within 71, the uniform 3x3 dyadic 7 = 7 plays the role 

of the metric with components (21)-(26), whereas the uniform 3 vector £ = T_ \ n plays the 
role of the vector with components (27)-(29). As a consquence of the uniform metric ap- 
proximation, henceforth we restrict our consideration to electromagnetic wavelengths which 
are small compared with the linear dimensions of the neighbourhood 1Z. 

Within the neighbourhood 1Z, the electromagnetic fields are decomposed via the three- 
dimensional Fourier transforms 



pOO poo poo 

E_(ct,r) — - III E(o;/c, k) exp [i(k • r — cut)] dm dk± dki 

C J —00 J~ 00 J —00 
poo poo poo 

H_{ct,r) — - III HO^ / c i k) exp [i(k • r — cut)] duj dk± dk 2 

C J —00 J —00 J —00 



(33) 



where the wavevector k is the Fourier variable corresponding to r, i — T and u is the 
usual temporal frequency. The wavevector component k 3 is determined by combining (33) 
with the Maxwell curl postulates (7) and then solving the resulting 4x4 eigenvalue problem 
[18]. We take k 3 G R and thereby consider only propagating (non-evanescent) plane waves. 

The associated rate of energy flow is calculated from the complex-valued phasors of the 
electric and magnetic fields, E(u>/c, k) and H(u>/c,k) respectively, by means of the time- 
averaged Poynting vector 1 

( P{cu/c,k) ) t = l - Re { E(u/ c, k) x H*(u/c,k) } . (34) 

If the phase velocity vector casts a negative projection onto the time-averaged Poynting 
vector then we have so-called negative phase velocity (NPV); i.e., NPV is signalled by 

k' <P(w/c,fc)) t < 0. (35) 

The key steps involved in deriving a general expression for the quantity k • (P(w/c,k) )t 
follow the corresponding steps for a simply moving medium which is isotropic dielectric- 
magnetic at rest [16], [19]. We begin by combining the Fourier decompositions (33) with the 
constitutive relations (8) and the Maxwell curl postulates (7). Thereby, we find 

px E(w/c, k) = ujfio^j • H(u/c, k) , (36) 
pxH(u/c,k) = — ue 7 • E(u/c, k) (37) 



^Re{-} denotes the real part while the asterisk denotes the complex conjugate. 
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where the vector 



p=k--L 

- c 



(38) 



has been introduced. Upon substituting (36) into (37) and carrying out some algebraic 
manipulations, we obtain the eigenvector equation 



|7| — V ' 7 * V 



L+PP'lf ' E (uj/c, k) = 



(39) 



and the corresponding dispersion relation 



p • 7 m p 



c J 



71 



0. 



(40) 



Here, / is the identity dyadic, and | 7 | denotes the determinant of 7. 
To solve (39), we exploit the orthogonality condition 

p • 7 • E(o;/c, k) = . 



(41) 



which follows from (40). The general solution to the electric phasor equation (39) is express- 
ible as the sum 



E(c<j/c, k) = A a {uj I c, k) e a (u / c, k) + A b (u / c, k) e b (u / c, k) 
where the unit vectors 



(42) 



7 • w 



17 m 

g 1 * (p x e a ) 
|7 _1 * (P x e a ) I j 



(43) 



are chosen in accordance with (41). The unit vector w is orthogonal to p, i.e., w % p — 0, 
but is otherwise arbitrary. For calculating the numerical results presented in the following 
section, we chose 



w = < 



2P 2 Z + ( PX +Pyf 
1 



(Pz,Pz,~Px -Py) for p 2 ^0 



(44) 



s=(Py,-Px,0) 



for p 2 = 



with p = (p x ,Py,Pz) in the Cartesian coordinate system. Initial and boundary conditions fix 
the complex-valued amplitudes A a ^{uj/c,k). The general solution to the magnetic phasor 
equation corresponding to (39) follows directly from (42) upon substitution into (36). 
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Finally, combining (39) and (36) yields the central result 

k'(P)t = (\ A a\\ 'X *e a + \A b \ 2 e b -g -e 6 ) k -g -p. 

For the Kerr metric, we observe that 

\g\=5 2 [R 4 + (a Thh z) 2 } 2 . 
Therefore, it follows that | 7 | > and the sufficient condition 



Pa<0 
Pb<0 



for NPV propagation arises, wherein 

Pa = (e a 'J m §a) (t '1 'P 

Pb= (e 6 '2 ' e 6 ) (k'X'P) 



(45) 



(46) 



(47) 



(48) 



The dispersion relation (40) yields two wavenumbers k = k^ for the arbitrarily oriented 
wavevector k = kk with k = (sin 6 cos 0, sin 6 sin 0, cos 6) . These are 



k «7 • f ± J [k «7 -f y -k «7 'k(t -7 -f - I7Q 



V 



k'^-k 



(49) 



As only propagating (non-evanescent) plane waves are considered here, we have k^ 1 G R. 
In order to establish whether k • (P) t < for the wavenumbers k^, we resort to numerical 
evaluations of (45) and (49). 



3 Numerical results 

The prospects for NPV propagation are investigated within the gravitational field of a ro- 
tating black hole, as characterized by the Kerr metric, by calculating P a ^ of (48) for all 
directions of wave propagation in a chosen 1Z. 

For compact presentation of our results, we considered surfaces of constant R for 



R+ ^ R 5; Rsx 



(50) 



The region represented by (50) is known as the ergosphere of the blackhole. In regions 
beyond the stationary limit, i.e., R > Rs + , the NPV conditions (47) were found not to be 
satisfied. This observation is consistent with the asymptotic behaviour of the Kerr metric: 
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in the limit R — > oo, Kerr spacetime converges to flat Minkowski spacetime which does not 
support NPV propagation in vacuum [19]. 

For the purposes of illustration, we chose a black hole with a rbh = m rbh a/3/4. We explored 
the feasibility of NPV propagation, in all directions denoted by k = (sin 9 cos 0, sin 9 sin 0, cos 9) , 
at (x,y,z) neighbourhoods on surfaces of constant R. In Figure 1, the surfaces 



(a) 


R = 


R + (i.e., the outer event horizon), 


(b) 


R = 


0.5R+ + 0.5R S+ , 


(c) 


R = 


0.25P+ + 0.75P S+ , 


(d) 


R = 


0.1R + + 0.9R S+ , and 


(e) 


R = 


0.05P+ + 0.95R S+ 



are mapped for the octant {x > 0, y > 0, z > 0}. The three-dimensional surface mesh is 
constructed from z heights calculated on a 90 x 90 grid of (x, y) points. At each (x, y, z) 
point, the wavevector k L = k ± k L , together with the parameters P a and P&, were computed at 
1° increments for 9 G [0°, 180°) and G [0°, 360°). If NPV was found to arise at a particular 
(9, <f>) orientation (as indicated by both P a < and P5 < 0) then the (St, y, z) grid point was 
given a 'score' of one unit." Thus, at each (x, y, z) grid point the maximum theoretical score 
is 180 x 360 x 2 = 129600. The NPV scores are represented by colours in Figure 1. We 
remark that the magnitude of u does not influence the signs of P a ^ (nor the NPV scores). 

It is clear from Figure 1 that NPV propagation is most prevalent at (x, y, z) neighbour- 
hoods close to the equator of the outer event horizon (i.e., at neighbourhoods in the vicinity 
of the plane z = 0). The density of NPV propagation decreases towards the x = y = pole 
of the R = R + region. We note that the ergosphere vanishes at the poles, as the outer event 
horizon and the stationary limit surface coincide. 

For ergospheric regions outside the outer event horizon, distributions of NPV propagation 
density are mapped in Figure l(b)-(e). On comparing the maps in Figure l(a)-(e), it 
becomes clear that, as the constant R surfaces approach the stationary limit surface, the 
densities of NPV propagation become progressively more concentrated around the z = 
equator. 

In order to investigate the influence of the angular velocity term a Ihh we repeated the 
calculations of Figure 1 but with a Ihh = m Ibh /3. The corresponding NPV propagation maps 
are presented in Figure 2, for the same surfaces of constant R as were used in Figure 1. As 
compared with a lhh = m Thh a/3/4, the outer event horizon and the stationary limit surface 
are closer together for a rbh = m rbh /3. The concentration of NPV propagation at the z = 
equator is evidently somewhat more pronounced in the case of the thinner ergosphere. It also 
noteworthy that the overall occurrence of NPV propagation is less common for the thinner 
ergosphere, as indicated by the generally lower NPV scores recorded in Figure 2 as compared 
with Figure 1. 

II If NPV occurs for both the k + and k~ wavenumbers for a particular (6, <j>) orientation then a score of 
two units was awarded. 
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From Figures 1 and 2 it is clear that the occurrence of NPV increases as the black 
hole rotation increases. It is therefore of interest to look at the scenario a rbh — + m rbh for 
which the occurrence of NPV is maximized. The NPV distribution maps corresponding to 
a rbh = 0.99m rbh are shown in Figure 3, for the same surfaces of constant R as were considered 
in Figures 1 and 2. In this maximized scenario, the regions which support NPV propagation 
are widely distributed throughout the ergosphere. However, we note the relative sparsity of 
NPV propagation in the vicinity of the black hole's axis of rotation. 

Extrapolating from Figures 1, 2 and 3, we infer that NPV propagation does not occur 

(i) in neighbourhoods along the black hole's axis of rotation, and 

(ii) in the limit as the black hole's rate of rotation becomes vanishingly small. 

We note that the outer event horizon and the stationary limit surface merge as the rotational 
rate approaches zero, and thereby the ergosphere disappears. 

On comparing the maximum theoretical NPV score of 129600 with the maximum scores 
in Figures 1, 2 and 3, it is apparent that the conditions for NPV propagation are satisfied 
for relatively few wavevector orientations in the ergosphere. To illustrate this point further, 
we considered a single representative neighbourhood within the ergosphere, say x = 0.2m rbh , 
y = 1.8m rbh and z = 0.3m Ibh , with a Thh = m Tbh a/3/4 (for which R = 1.63 m lhh ). We 
computed P ab for all orientations of k. At orientations where P a ^ < a score of one unit 
was recorded. The scores are plotted against the spherical polar coordinates 9 and of k in 
Figure 4. Orientations at which P ah < arises from the k + wavenumber are coloured red, 
whereas those orientations at which P ab < occurs for the k~ wavenumber are coloured 
green. The k + NPV regions and the k~ NPV regions are observed in Figure 3 to occupy a 
relatively small portion of the (9, (f>) plane. 

4 Concluding remarks 

By using the Kerr metric in the Tamm formulation of electromagnetics in gravitationally af- 
fected vacuum, we have mapped the occurrence of negative-phase velocity (NPV) planewave 
propagation in the ergosphere of a rotating black hole. Our numerical results indicate that 

• NPV propagation does not appear possible outside the stationary limit surface; 

• rotation of a black hole is required for NPV propagation; 

• NPV propagation is impossible on the axis of rotation of the black hole; 

• regions of NPV propagation are concentrated at the equator of the ergosphere, and are 
less common towards the polar regions; and 

• an increase in angular velocity exaggerates the NPV concentration at the equator. 
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In particular, the importance of the rotational character of the black hole metric should be 
stressed. In the absence of rotation the black hole metric does not support NPV propagation. 
The infeasibility of NPV propagation for (i) a non-rotating black hole; and (ii) along the axis 
of rotation of a black hole, is consistent with our earlier analytic results [5, 6]. In connection 
with this issue, we note that NPV propagation can however arise in Schwarzschild-de Sitter 
spacetime [20]. 

Let us comment upon the applicability of the NPV condition (47). Suppose that the 
neighbourhood 1Z has representative linear dimensions given by p. In order for the nonuni- 
form metric g a p to be approximated by the uniform metric g a p throughout 71, it is necessary 
that p be small compared to the radius of curvature of the g a p spacetime. A convenient 
measure of the inverse radius of spacetime curvature squared is provided by the nonzero 
components of the Riemann tensor. These components are of the order of R~ 2 for the Kerr 
metric [21]. Therefore, we have p <C R. Since the neighbourhood TZ is also required to be 
large compared with electromagnetic wavelengths, as given by 2ir/\k\, we see that (47) holds 
in the regime 

(51) 

For black holes of the mass of our sun, the condition (51) translates to wavelengths of the 
order of 100 metres or less, whereas for supermassive black holes of the type believed to lie 
at the centre of galaxies the NPV condition (47) holds for wavelengths of the order of 10 7 
kilometres or less, typically [22]. 

The NPV phenomenon reported in the preceding sections concerns the propagation of 
electromagnetic plane waves in the ergosphere. The issue of possible generalizations to scalar 
and/or gravitational waves is a matter for future study. However, we note in this context 
that the Teukolsky equation yields general (nonplanar) solutions which support negative 
wavenumbers for scalar, electromagnetic and gravitational waves in the ergosphere of a Kerr 
black hole [23]-[25]. 

It is of interest to speculate whether the NPV propagation demonstrated here may arise 
for so-called acoustic black holes [26], typically studied in the cylindrical geometry [27]. 
While cylindrical acoustic black holes can support superradiant scattering analogously to 
the equatorial slice of the Kerr black hole [27, 28] , the propagation of NPV acoustic waves 
for acoustic black holes has yet to be established. 

We close by emphasizing that the description of NPV presented herein arises from the 
classical treatment of electromagnetic waves in curved spacetime, within the short wavelength 
regime indicated by (51). This NPV phenomenon is distinct from superradiant scattering in 
the ergosphere at long wavelengths 
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Figure 1: NPV maps for a rbh = m rbh a/3/4. See the text for an explanation of the colour 
mapping, (a) At the surface R = R + (i.e., at the outer event horizon) (top left); (b) At 
the surface R = 0.5R+ + 0.5R S+ (top right); (c) At the surface R = 0.25R+ + 0.75R S+ 
(middle left); (d) At the surface R = 0.1R + + 0.9Rs + (middle right); and (e) At the surface 
R = 0.05R+ + 0.95R S+ (bottom) . 
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Figure 2: NPV maps for a Ihh = m Ihh /3. See the text for an explanation of the colour 
mapping, (a) At the surface R = R + (i.e., at the outer event horizon) (top left); (b) At 
the surface R = 0.5i?+ + 0.5R S+ (top right); (c) At the surface R = 0.25R+ + 0.75R S+ 
(middle left); (d) At the surface R = 0.1-R + + 0.9Rs + (middle right); and (e) At the surface 
R = 0.05R+ + 0.95R S+ (bottom). 
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Figure 3: NPV maps for 

<^rbh — 0.99T7i rbh . See the text for an explanation of the colour 
mapping, (a) At the surface R = R + (i.e., at the outer event horizon) (top left); (b) At 
the surface R = 0.5R+ + 0.5R S+ (top right); (c) At the surface R = 0.25R+ + 0.75i2 s+ 
(middle left); (d) At the surface R = 0.1R + + 0.9Rs + (middle right); and (e) At the surface 
R = 0.05R+ + 0.95R S+ (bottom). 
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Figure 4: NPV scores plotted as heights, against the k orientational angles 9 and 0, in the 
neighbourhood of (x = 0.2m rbh , y = 1.8m rbll , z = 0.3m lbh ), with a rbh = 77T. rbh ^3/4. NPV 
arising from the k + wavenumber is coded red; NPV arising from the k~ wavenumber is coded 
green. 
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